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dynamics 
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New generalized equations of motion for the Weber-Clebsch potentials that describe both the 
Navier-Stokes and MHD dynamics are derived. These depend on a new parameter, which has 
dimensions of time for Navier-Stokes and inverse velocity for MHD. Direct numerical simulations 
(DNS) are performed. For Navier-Stokes, the generalized formalism captures the intense recon- 
nection of vortices of the Boratav, Pelz and Zabusky flow, in agreement with the previous study 
by Ohkitani and Constantin. For MHD, the new formalism is used to detect magnetic reconnec- 
tion in several flows: the 3D Arnold, Beltrami and Childress (ABC) flow and the (2D and 3D) 
Orszag-Tang vortex. It is concluded that periods of intense activity in the magnetic enstrophy are 
correlated with periods of increasingly frequent resettings. Finally, the positive correlation between 
the sharpness of the increase in resetting frequency and the spatial localization of the reconnection 
region is discussed. 

PACS numbers: 47.10.-g, 47.ll.-j, 47.32. C-, 47.65.-d 



I. INTRODUCTION 

The Eulerian-Lagrangian formulation of the (inviscid) 
Euler dynamics in terms of advected Weber-Clebsch po- 
tentials [l| was extended by Constantin [2[ to cover the 
(viscous) Navier-Stokes dynamics. Ohkitani and Con- 
stantin (OC) 3] then performed numerical studies of this 
formulation of the Navier-Stokes equations. They con- 
cluded that the diffusive Lagrangian map becomes non- 
invertible under time evolution and requires resetting for 
its calculation. They proposed the observed sharp in- 
crease of the frequency of resettings as a new diagnostic 
of vortex reconnection. 

We were able to recently complement these results, us- 
ing an approach that is based on a generalized set of equa- 
tions of motion for the Weber-Clebsch potentials that 
turned out to depend on a parameter r which has the 
unit of time for the Navier-Stokes case Q (the MHD case 
is different, sec below Section III A 3p . The OC formula- 
tion is the (singular) r — » limit case of our generalized 
formulation. Using direct numerical simulations (DNS) 
of the viscous Taylor-Green vortex |5J we found that for 
t/0 the Navier-Stokes dynamics was well reproduced at 
small enough Reynolds numbers without resetting. How- 
ever, performing resettings allowed computation at much 
higher Reynolds number. 

The aim of the present article is to extend these results 
to different flows, both in the Navier-Stokes case and in 
magnetohydrodynamics, and thereby obtain a new diag- 
nostic for magnetic reconnection. Our main conclusion is 
that intense reconnection of magnetic field lines is indeed 
captured in our new generalized formulation as a sharp 
increase of the frequency of resettings. Here follows a 



summary of our principal results. 

We first derive new generalized equations of motion 
for the Weber-Clebsch potentials that describe both the 
Navier-Stokes and MHD dynamics. Performing DNS of 
the Boratav, Pelz and Zabusky flow [6|], that was pre- 
viously used by Ohkitani and Constantin Q, we first 
check that our generalized formalism captures the in- 
tense Navier-Stokes vortex reconnection of this flow. We 
demonstrate the reconnection of vortices is actually oc- 
curring at the instant of intense activity in the enstro- 
phy, near the lows of the determinant that trigger the 
resettings. We then study the correlation of magnetic 
reconnection with increase of resetting frequency by per- 
forming DNS of several prototypical MHD flows: the 3D 
Arnold, Beltrami and Childress (ABC) flow and the 
Orszag-Tang vortex in 2D [8[ and 3D Q . 



II. THEORETICAL FRAMEWORK 

A. General Setting 

1. Weber-Clebsch representation for a class of evolution 
equations 

Let us consider a 3D vector field Z depending 
on time and (3-dimcnsional) space, with coordinates 
(x ,x , x 3 , t). Assume Z satisfies an evolution equation 
of the kind: 



L»Z 
~~Dt 



-X7P+J2u a VZ a + nAZ (1) 



V-Z = 0, 



(2) 



where greek indices a, (3 denote vector field components 
running from f to 3, u is a given 3D velocity field and 
we have used the convective derivative defined by 
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In the following sections, two different cases will be 
considered. In section Hi Bl (Navier-Stokes case) the vec- 
tor field Z will correspond to the velocity field u, whereas 
in section lll CI (MHD case) it will correspond to the mag- 
netic vector potential A. 

Let us first recall that performing a change from La- 
grangian to Eulerian coordinates on the Weber transfor- 
mation [101 ] leads to a description of the Euler equations 
as a system of three coupled active vector equations in a 
form that generalizes the Clebsch variable representation 

Q. 

Our starting point will be to apply this classical Weber- 
Clebsch representation to the field Z: 



Z = ^A^V/Z-V^, 



(3) 



i=l 



where each element of the 3 pairs of Weber-Clebsch po- 
tentials (A 1 ,//), i = 1,2, 3 is a scalar function. 

Performing a variation on the Weber-Clebsch represen- 
tation ([3]) yields the relation 
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5Z = V (<5A l V// - c5//VA ! ) - V 



=i 



3 

E 

(=1 



<VA 4 



( 4 ) 

where the symbol 5 stands for any (spatial or tempo- 
ral) partial derivative. Taking into account the identity 
[V, -gj] = (Vu) ■ V, it is straightforward to derive from 
(|U) the following explicit expression for the convective 
derivative of the vector field Z: 



J 



- W^v-^w,- 



Dt 
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a=l 

r 



(5) 



2. Equations of motion for the potentials 

Following steps that are similar to those presented in 
our previous paper [J] , we now derive a system of equa- 
tions of motion for the Weber-Clebsch potentials ([3|) that 
is equivalent to the original equation |1| . If we use the 
RHS of equation (p} to replace the LHS of our general 
identity (O, the resulting relation can be solved for the 
time derivative of the potentials: 



DX 
~Dt 

~~Dt 



= kAX 1 + L l [X, fi] 
= kA^ + M l [X,fj] 



Here U,M % obey the linear equation 



E (#v 



M l VA' 



VG. 



(6) 
(7) 

(8) 



The "divergence-less gauge" ([2]) allows one to express 
in terms of X 1 and fi 1 , as the solution of the linear 
equation 



A^ = ]TV-(A*V^)- 



(11) 



Thus there is no need to solve equation (|T0|) for the field 
0, since this equation is identically satisfied when <f> is 
determined by eq. (fTTj) . 

Equation 1(5)1 above is a system of 3 linear equations for 
the 6 unknowns L 1 , M % . When n = there is a simple 
solution to ([5]): L l = M l = G = 0. In this case the evo- 
lution equations ([6]) and ([7]) represent simple advection. 



Moore-Penrose solution and minimum norm 



where 



f = 2K EE^ AI5 « v ^ 



i=X a=l 



(9) 



and G[A,/i] is an arbitrary scalar related to the non- 
unique separation of a gradient part in eq. (O : 



The linear system ((8]) is underdetermined (3 equations 
for 6 unknowns). In order to find a solution to the sys- 
tem we need to impose extra conditions. Since L l ,M l 
appear in the equations on an equal footing, it is natural 
to supplement the system by a requirement of minimum 
norm, namely that 



td±-P = Y \*M* -G-u-Z. 



(10) 



J2(L l L* 



T -' 1 M l M i ') 



(12) 



be the smallest possible (this is the so-called general 
Moore-Penrose approach [ll], Qj], [l3| , see also our previ- 
ous paper J4|]). The parameter r has physical units equal 
to [M/L]. Using eqs.([H]),([7]) these are the units of [///A]. 
It will turn out (see equation (|20[) below) that [/j] = L 
(length) and this implies from eq.© that [A] = [Z]. 
Therefore the units of r are 

L 



In the Navier-Stokes case (section lriBJ) [Z] = [u] = LT^ 1 
and thus [r] — T, whereas in the MHD case (section lllCp 
[Z] = [A] = I?T- X and thus [r] = TL~ X . 

The Moore-Penrose solution to ([5]) , that minimizes the 
norm (J12p. is explicitly given in equations (A6,A7) of 
reference [J]. Inserting this solution in ©,(17]) we finally 
obtain the explicit evolution equations 



DX 1 
~Dt 
Drf 

~Dt 



kAX 1 + V// ■ 
kA// - t 2 VA' 



f-VG 



f-VG 



(13) 
(14) 



where f is given in eq.©, the dot product denotes ma- 
trix or vector multiplication of 3-dimensional tensors, and 
HU 1 is the inverse of the square symmetric 3x3 matrix 
H, defined by its components: 



H. 



a/3 



^2(r i a a x i dfiX i + d a n i d l 



w 



(15) 



i=l 



These evolution equations together with the particular 
choice for the arbitrary function G (see equation (All) 
of reference |4|]) 



G = A _1 V-f, 



(16) 



is our new algorithm. 

In the Navier-Stokes case, we showed in a previous pa- 
per [J that the limit r — > corresponds to the approach 
used by Ohkitani and Constantin Q . In the general case 
(Navier-Stokes as well as MHD), we remark that the ma- 
trix H (see equation (fT5|) ) can be written (using obvious 
notation) as H = (V/x) • (V/x) T + r 2 (VA) • (VA) T , which 
has a very simple structure in the limit r — > . Be- 
cause the condition det(V/i) = is generically obtained 
at lower codimension than the condition det H = 0, the 
limit r — > is singular. 



B. Navier-Stokes equations 

The standard incompressible NS equations can be writ- 
ten in the form: 

— j = -V U+-|u| 2 J +^u tt Vu a +^Au 

V-u = 0, 

which is indeed of the general form |T]), (2J with Z = u, 
k = v and P = p+ jluj 2 . 



C. MHD equations 

The standard incompressible MHD equations for the 
fluid velocity u and the induction field b, expressed in 
Alfvenic velocity units, can be written in the form: 

— -^ = -Vp + i/Au+(Vxb)xb (17) 



— = (b-V)u + 7?Ab 

V-u = 
Vb = 0, 



(18) 
(19) 



where v and r\ are the viscosity and magnetic resistivity, 
respectively. 

We introduce the vector potential in the Coulomb 
gauge: 

b = V x A 

V-A = 0. 

Using the identity V x (J2 a =i u a^ A a — (u • V)A) = 
(b • V)u — (u • V)b — (V • u)b and the incompressibility 
condition (fi"9| , eq. (fTS"]) can be written as 

— = -Vc+^u Q V/l«+7/AA, 

which is indeed of the general form {T]) with Z = A, k = r\ 
and P = c. 



III. NUMERICAL RESULTS 
A. Implementation 

1. Initial conditions in pseudo- spectral method 

Spatially periodic fields can be generated from the 
Weber-Clebsch representation $5§ by setting 

// = ^ + M' P , (20) 

and assuming that fi 1 and the other fields A 1 and (f> ap- 
pearing in ([3]) are periodic. Indeed, any given periodic 
field Z can be represented in this way by setting 



/'p 
A' 



= 

= z l 

= 0. 



(21) 
(22) 
(23) 



Note that the time independent non-periodic part of fi l 
of the form given in (|20p is such that the gradients of 
H l are periodic. It is easy to check that this representa- 
tion is consistent with the generalized equations of mo- 
tions (|13ll4p . We chose to use standard Fourier pseudo- 
spectral methods, both for their precision and for their 
ease of implementation [14| . 



2. Resettings and reconnection 

Following Ohkitani and Constantin [3], we now de- 
fine resettings. Equations ^T\i , (HU) and (f!?o) are used 
not only to initialize the Weber-Clebsch potentials at the 
start of the calculation but also to reset them to the cur- 
rent value of the field Z, obtained from © and (fTTj) . 
whenever the minimum of the determinant of the matrix 
(fT5} falls below a given threshold 

detH<e 2 . 

It is possible to capture reconnection events using re- 
settings. The rationale for this approach is that recon- 
nection events are associated to localized, intense and in- 
creasingly fast activity which will drive the potentials to 
a (unphysical) singularity in a finite time. One way to de- 
tect this singularity is via the alignment of the gradients 
of the potentials, which leads to the vanishing of det HI 
at the point (s) where this intense activity or 'anomalous 
diffusion' is taking place. Now, the time scale of this 
singularity is much smaller than the time scale of the re- 
connection process itself Q, so when detH goes below 
the given threshold and a resetting of the potentials is 
performed, the anomalous diffusion starts taking place 
again, more intensely as we approach the fastest recon- 
nection period, driving the new (reset) potentials to a 
new finite-time singularity, in a time scale that decreases 
as we approach this period. Therefore, successive reset- 
tings will be more and more frequent near the period of 
fastest reconnection, and that is what we observe in the 
numerical simulations. This procedure will be used to 
capture reconnection events in particular flows in both 
the Navier-Stokes case (Z = u, Section IIIIBI) and the 
MHD case (Z = A, Sections HHII! and |TTlC3| . 



B. Navier-Stokes case: BPZ Flow, resettings and 
reconnection 



1. Comparison of Weber-Clebsch algorithm with DNS of 
Navier-Stokes 

In order to characterize the precision of the r ^ 
Weber-Clebsch algorithm, we now compare the veloc- 
ity field Z = u obtained from §S§ and (fTTj) , by evolving 
the Weber-Clebsch potentials using (|T5j) -" ([TB)l . with the 
velocity field obtained independently by direct Navier- 
Stokes evolution from the BPZ initial data. 

More precisely, we compare the associated kinetic en- 
strophy £l(t) = X^j. fc 2 i£(A;,i) where the kinetic energy 
spectrum E(k, t) is defined by averaging the Fourier 
transform u(k',t) of the velocity field $5§ on spherical 
shells of width Afc = 1, 



E{k,t) = \ 



E 



|u(k',t) 



fe-Afe/2<|k'|<fe+A/c/2 




Ohkitani and Constantin (OC) [3[ used a flow that ini- 
tially consists of two orthogonally placed vortex tubes 
that was previously introduced in Boratav, Pelz and 
Zabusky (BPZ) [6] to study in detail vortex reconnec- 
tion. Our previous numerical study of the generalized 
Weber-Clebsch description of Navier-Stokes dynamics [4j 
was performed using the Taylor-Green vortex, a flow in 
which vorticity layers are formed in the early stage, fol- 
lowed by their rolling-up by Kelvin-Helmholtz instability 
[15| . It can be argued [3} that cut-and-connect type re- 
connections are much more pronounced in the BPZ flow 
than in the Taylor-Green flow. In this section we present 
comparisons, performed on the BPZ flow, of our t =/= 
generalized algorithm with direct Navier-Stokes simula- 
tions and with OC original approach. The potentials 
are integrated with resettings in resolution 128 3 for a 
Reynolds number of R = 1044, which is the one used by 
BPZ and OC. 

The BPZ initial data is explicitly given in [6|]. 



FIG. 1: Navier-Stokes case: BPZ Flow. Temporal evolution 
of kinetic enstrophy fl for a Reynolds number of R — 1044 
with t — 0, 0.01 and 0.1 (+, o and x). The solid line comes 
from a direct numerical simulation (DNS) at resolution 128 3 . 



Figure Q] shows that the kinetic enstrophy is well re- 
solved, independently of the choice of the parameter r. 



2. Time between resettings as a method for reconnection 
capture 

In this section we study the influence of the parame- 
ter r on the temporal distribution of the intervals Atj = 
tj —tj—i between resetting times tj, at fixed value of the 
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FIG. 3: Visualization of vorticity lu for the BPZ flow (same 
conditions than in Figures [l] and [2}. Note the change of topol- 
ogy of the vortex tubes before (left, t = 3.2, uVax = 20) and 
after (right, t — 7.1, cj max = 15) the reconnection process. 
Isosurfaces colors: orange: 6, yellow: 9, green: 12 and blue: 
16 (color online). 



FIG. 2: Navier-Stokes case: BPZ flow. Temporal evolution of 
resetting interval At for r = 0, 0.01 and 0.1 (o, □ and +), the 
triangles correspond to the simulation performed by Ohkitani 
and Constantin. 



resetting threshold e 2 = 0.1. Using the same Reynolds 
number and resolution that was used to create Fig. [1] 
Figure [5] is a plot of At as a function of time, for simu- 
lations with different values of r. In the same figure we 
also show the corresponding At for a replica of the sim- 
ulation performed by OC, that is in excellent agreement 
with our general case. 

We see that, independently of r, there are sharp min- 
ima in At during the periods of maximum enstrophy (see 
Fig. [1} . Inspection of figure [3] demonstrates that the 
deepest minimum corresponds in fact to the time when 
reconnection is taking place. The main tubes in the left 
and right figures are isosurfaces of vorticity correspond- 
ing to 60% of the maximum vorticity, which is attained 
inside each of the main tubes. 

Figured] (left) shows that the spatial region where the 
determinant det EI goes below the threshold before each 
resetting corresponds to a small, localized neighborhood 
between the main interacting vortices. This region is seen 
in the right figure as a bridge connecting the two vortices: 
this bridge is an isosurface of vorticity corresponding to 
73% of the maximum vorticity, which is attained inside 
the bridge. The main tubes correspond to isosurfaces of 
30% of the maximum vorticity. Note that this behavior 
of the determinant det HI is also true for any value of 
t (data not shown), confirming in this way the original 
rationale for the study of reconnection with the aid of 
resettings. 

Figures [3] and [J] were made using the VAPOR [la 
fl3l visualization software. 



C. MHD Flows 

In this section we study MHD flows with simple initial 
conditions. The magnetic potential Z = A is obtained 





FIG. 4: Visualization of the determinant of the matrix (|15[) 
(left) and vorticity (right, with the same color map than in 
Figure [3]) at reconnection time t — 4.7, w max = 43.6 (see 
Figure [3]), for r = 0.01. The region where the determinant 
triggers resetting is within the displayed blue isosurface at 9 
times the triggering level e 2 (color online). 



in terms of the Weber-Clebsch potentials from ([3J and 
(fTTj) . and the Weber-Clebsch potentials are evolved using 
equations ([I3" ]) -([T6" ]) . 

We treat the evolution of the velocity field in two differ- 
ent ways: (i) As a kinematic dynamo (ABC flow, Section 
IIIIC ip . where the velocity is kept constant in time; (ii) 
Using the full MHD equations (Orszag-Tang 2D and 3D, 
Sections IIIIC 21 and IIIIC 31) , where the velocity field is 
evolved using the momentum equation (jTTJ) . 

To compare with DNS of the induction equation (JT8J) 
for the magnetic field we proceed analogously as in the 
Navier-Stokes case. We compare the magnetic enstrophy 
[18| Q m (t) = J2k k 2 E m (k,t), where the magnetic energy 
spectrum E m (k,t) is defined by averaging the Fourier 
transform b(k', t) of the magnetic field b = VxA (with 
A given by ([3])) on spherical shells of width Afc = 1, 



E m (k,t) 



fc-A/c/2<|k'|<fc+Afc/2 



b(k',t)| 



Note that magnetic dissipation is the square current. 
Resettings will be performed with a resetting threshold 



e 2 = 0.1. We have checked that e 2 = 0.4 and e 2 = 0.025 
give results that vary only slightly (figures not shown). 
This is an evidence of the robustness of the resetting 
method and a validation of the rationale for the use of 
resettings to diagnose reconnection. 



1. Kinematic dynamo: ABC Flow 



We have used the ABC [7] velocity: 



u x = Bq cos koy + Co sin koz 
u y = C'o cos koz + Ao sin kox 
u z = Aq cos kox + Bq sin koy , 



with k — 2 and A = B a 
magnetic seed that reads 



C, 







1. We used an initial 



A a 

A/ 

.4, 






do sin x sin y . 



The magnetic resistivity has been chosen as r/ = 1/12 
and we have set do — 1/100 for simplicity (its value is 
unimportant in the kinematic dynamo). 

Runs with resettings are compared for different values 
of the parameter r. It is seen in Fig. Othat the magnetic 
enstrophy Q, m is well resolved for each case, at resolution 
128 3 . 

The resettings are quite regular in time and indeed 
they slow down as time goes by, at a regular rate which 
decreases with increasing resolution (figure not shown). 
There is no increase in the resetting frequency. This be- 
havior is consistent with the monotonic behavior of the 
magnetic enstrophy and with the absence of localized or 
intense activity of the magnetic field. 




FIG. 5: Temporal evolution of magnetic enstrophy Q. m for the 
ABC flow (kinematic dynamo) with r = and 1 (o and +), 
with a resolution of 128 3 and the constants v — and r\ — 
1/12. The solid line comes from a direct numerical simulation 
(DNS) of the induction equation for the magnetic field. 



2. Full MHD equations: 2D Orszag-Tang Vortex 

In the rest of the paper, the full MHD equations of mo- 
tion are integrated. The momentum equation for the ve- 
locity (fTTj) is integrated together with the Weber Clebsch 
evolution equations (fT"3)) - (fl"6|) where the magnetic poten- 
tial Z = A is obtained from $5§ and ([Til .. 

We have chosen the following initial data for the 2D 
Orszag-Tang (hereafter, OT) vortex [g|: 



Ux 


= 


—2smy 


Uy 


= 


2 sin a; 


u z 


= 


o, 


A x 


= 





Ay 


= 





A z 


= 


2 cos x cos 2y 



The OT vortex has a magnetic hyperbolic X-point lo- 
cated at a stagnation point of the velocity, and is a stan- 
dard test of magnetic reconnectiom both in two dimen- 
sions [l9| and in three dimensions [20]! see below section 
IIIIC 31 

We compare runs with resettings for different values 
of the parameter r. Figure [6] shows that the magnetic 
enstrophy is well resolved in resolution 128 2 . 

Figure [7] shows the time between resettings as a func- 
tion of time, for runs performed with different values of 
t. It is apparent from the figure that there are periods 
of frequent resettings which coincide with the periods of 
high magnetic enstrophy from Fig. [5] This is a robust 
evidence of the utility of the resetting approach for 2D 
magnetic reconnection. 

We have also simulated the Orszag-Tang vortex in the 
so-called 2.5D setting [21[ (see also the DiPerna-Majda's 
construction [22|]), defined by the same initial data as the 
above 2D Orszag-Tang vortex, but with A x = sin y and 
A y = — sinx. We obtained (data not shown) a behavior 
of the resetting frequency which was very similar to that 
of the 2D case. 



3. Full MHD equations: 3D Orszag-Tang Vortex 

For the 3D Orszag-Tang vortex 0] the initial magnetic 
potential reads 

A x = Co (cosy — cos z) 
Ay = Co (—cos x + cos z) 
A z = Co (cosx + cos2y) , 

with Co = 0.8. The initial velocity is given by 

u x — — sin y 



u z = 0. 

As in the 2D case, we compare runs with resettings for 
different values of the parameter r. Figure [5] shows 




FIG. 6: Temporal evolution of magnetic enstrophy Q m for 
Orszag-Tang in 2D for t — 0, 0.01 and 1 (o, □ and o) with 
a resolution of 128 2 and r\ = v — 0.005. Solid line: Direct 
numerical simulation of MHD equations. 
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< 1 



0.5 




FIG. 7: Temporal evolution of At for r = 0, 0.01, 0.1 and 1 
(o, □, o and A), for a simulation of Orszag-Tang in 2D with 
v — r\ = 0.005 and a resolution of 128 2 . 



that the magnetic enstrophy is well resolved in resolu- 
tion 128 3 , and Fig. [9] shows the time between resettings 
as a function of time. Again the periods of frequent re- 
settings coincide with the periods of high magnetic en- 
strophy from Fig. [8j proving the utility of the resetting 
approach for 3D magnetic reconnection. 



IV. CONCLUSIONS 

We have shown that the generalized Weber-Clebsch 
evolution equations allow to study reconnection events 
for both Navier-Stokes and MHD dynamics. We have 
checked for the Navier-Stokes BPZ flow that reconnec- 
tion events can be viewed as periods of fast and localized 
changes in the geometry of the Weber-Clebsch potentials, 
leading to more and more frequent resetting of the po- 



G20- 




FIG. 8: Temporal evolution of magnetic enstrophy Q m for 
Orszag-Tang in 3D for r = 0, 0.1 and 1 (o, □ and o) with 
a resolution of 128 3 and r\ = v = 0.005. Solid line: Direct 
numerical simulation of MHD equations. 
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FIG. 9: Temporal evolution of At for r = 0, 0.1 and 1 (o, □ 
and o), for a simulation of Orszag-Tang in 3D. 



tentials. 

We have applied the new generalized Weber-Clebsch 
evolution equations to the study of magnetic reconnec- 
tion in MHD. Taking as examples both the 2D and 3D 
Orszag-Tang vortices, we show a correlation of the re- 
connection events (associated to periods of high magnetic 
dissipation) with the periods of fast changes in the geom- 
etry of the Weber-Clebsch potentials, leading to frequent 
resettings of the potentials. 

However, unlike the case of BPZ reconnection, in this 
case the frequency of resettings does not have a sharp 
peak but a smeared one. Notice that, in the Navier- 
Stokes case, the corresponding frequency of resettings 
for the Taylor-Green vortex has also a mild peak. [J| 
One can argue that the 2D and 3D Orszag-Tang flows 



are more similar to Taylor-Green than to BPZ. Indeed, 
both Orszag-Tang and Taylor-Green have initial condi- 
tions with just a few Fourier modes, therefore they are 
extended spatially, whereas the BPZ initial condition is 
spatially localized (two orthogonal vortex tubes). 

This wide spatial extent of the vorticity in both 
Orszag-Tang and Taylor-Green vortices, as opposed to 
the localized extent of BPZ, might be the reason for the 
mildness in the shape of the minimum of the time be- 
tween resettings. In both spatially extended cases one 
expects reconnection events to happen in relatively dis- 
tant places at similar times, as opposed to the BPZ very 
localized cut-and-connect type of reconnection. In terms 
of the singularities of the Weber-Clebsch potentials and 
associated resetting, we should observe (to be studied in 
detail in future work) that the set of points where det H 
goes below the threshold consists of an extended region, 



as opposed to BPZ where we have confirmed that these 
points belong to a very localized region in space. Con- 
sequently, the widely distributed events that lead to re- 
setting in Orszag-Tang and Taylor-Green configurations 
would tend to be less correlated in time, leading to the 
smearing of the minimum of the curve for the time be- 
tween resettings, which would otherwise be very sharp if 
the events were more localized and therefore more corre- 
lated in time. 
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